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Abstract 



This work focuses on Exchangeable Occupancy Models (EOM) and their relations with 
the Uniform Order Statistics Property (UOSP) for point processes in discrete time. As our 
main purpose, we show how definitions and results presented in Shaked, Spizzichino and 
Suter (2004) can be unified and generalized in the frame of occupancy models. We first 



show some general facts about EOM's. Then we introduce a class of EOM's, called Ai^- 
modcls, that is used for generalizing the notion of Uniform Order Statistics Property. For 
processes with this property, we prove a general characterization result. Our interest is 
also focused on properties of closure w.r.t. some natural transformations of EOM's. 
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1 Introduction 

The so called occupancy distributions give rise, as well known, to a class of multivariate models 
useful in the description of randomized phenomena. The name "occupancy" comes from the 
interpretation in terms of particles that are randomly distributed among several cells. In 
particular, three classical examples, related to as many well-known physical systems, belong 



to this class: Maxwell-Boltzmann, Bose-Einstein and Fermi-Dirac model, see Feller (1968). In 



the Maxwell-Boltzmann (MB) statistics the capacity of each cell is unlimited, and the particles 
are distinguishable. In the Bose-Einstein (BE) statistics, the capacity of each cell is unlimited 
but the particles are indistinguishable. In the Fermi Dirac (FD) statistics, the particles are 
indistinguishable, but cells can only hold a maximum of one particle. All these three statistics 
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assume that the cells are distinguishable. These models are attractive for many reasons. First 
of all they have wide range of applications in Sciences, Engineering and also in Statistics, as 



pointed out from Charalambides (2005, Chapters 4 and 5) and Gadrich and Ravid (2011). 



Moreover Mahmoud (2009) provided an interpretation of occupancy distributions in terms 



of Polya Urns, which are very flexible and applicable to problems arising in various areas; 



e.g. Clinical Trials (see Crimaldi and Leisen (2008) for some references), Economics (see 



Aruka (2011)) and Computer Science (see Shah, Kothari, Jayadeva and Chandra (2011)). 
From a probabilistic and combinatoric point of view, the three fundamental models (MB, 
BE and BS) have many interesting properties. Indeed they are, in particular, exchangeable 
and this is a basic remark for our aims. This paper is in fact concentrated on the theme of 
Exchangeable Occupancy Models (EOM) and their relations with the Uniform Order Statistics 
Property (UOSP) of counting processes in discrete time. As one main purpose of ours, we 



show that some notions and results given in Shaked, Spizzichino and Suter (2004 2008) 



admit completely natural generalizations in the frame of EOM's. 

We start proving some general properties of the EOM's. In particular, we show that the ex- 
changeability property is maintained under remarkable types of transformations of occupancy 
models. Then, we introduce the notion of A4^ a )-models, a relevant sub-class of EOM's that 
turns out to have an important role in our derivations. In the final part of the paper, we 
study discrete-time counting processes satisfying the generalized UOSP; in other words, for 
which the joint distribution of the jump amounts is distributed according to a .M^-model, 
conditionally on a fixed number of arrivals up to time t. In particular, for these processes, 
several characterizations are proved. 

More in detail, the outline of the paper is as follows. Our central results will be presented in 
the Section [5j where we show natural extensions of the results proved in Shaked, Spizzichino 
and Suter | ( |2004| [2008] ) . Section [2] is devoted to some preliminary arguments about occupancy 
models and Section [3] presents some specific aspects of the EOM's. In Section [4] we define the 
class of occupancy -models and study some related properties. 



2 Occupancy models and order statistics of discrete variables 

For fixed n = 2, 3, . . . and r = 1, 2, . . . , let A nr be the set defined by 



x = (xi, . . . , x n ) : Xj = 0, 1, . . . , r and Xj 

i=i 



As it is well-known (see e.g. Feller (1968)) the cardinality of A n ^ r is 



| An,', 



n + r — 1 
n — 1 
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Starting from the classical scheme of r particles that are distributed stochastically into n 
cells, we consider the random vector X = (X±, X2, . . . , X n ) where, for j = 1, 2, . . . , n, Xj is 
the {0, 1, . . . , r}-valued random variable that counts the number of particles fallen in j-ih cell. 
Xi, X2, . . . , X n are called occupancy numbers and the joint distribution of (X±, X2, . ■ . , X n ), 
describing the probabilistic mechanism of assignment of the particles to the cells, is called 
an occupancy model. Since the vector (X\, X2, . ■ ■ , X n ) takes its values in the set A Utr , an 
occupancy model is then a probability distribution on A n ^ r . 

Let now B r ^ n denote the set 

B r ,n '■= {u = (u±, U2, • • • , u r ) : Ui = 1, 2, . . . , n and u\ < U2 < • • • < u r } 
and consider the mapping (p : B r ^ n — > A n ^ r defined as 

ip(u) = (ipi(u),ip 2 {u), v?n(u)) , 

with 



<M U ) = ^2 1 {ui=j} > for j = 1,2,..., n. 

i=i 

This is a one-to-one correspondence and then, for the cardinality of B rjJl , we have 
As to ip = tp~ l : A n: , r — > B r:n , we can write 

V>(x) = (V'l(x),^ 2 (x),...,Vr(x)), 



n + r — 1 
n- 1 



with 

Y>i(x) = min I s 



3=1 



> , for i = 1, 2, . . . , r. 

We can thus consider the random vector U = (U\, . . . , U r ) defined by 

U = ^(X). (1) 

For (ui, . . . , u r ) £ B rjn , one then has 

P{Ui =U 1 ,...,U r =U r } = P{Xi = (fi(u), ...,X n = <^n(u)}. 

Let V (A ntr ) and V (B r ^ n ) respectively denote the family of probability distributions on A n>r 
and the family of probability distributions on B r ^ n . In view of the above one-to-one corre- 
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spondence between A n>r and B r<n , we can consider the induced (one-to-one) correspondence 
between V (A n ^ r ) and V (B r ^ n ). More precisely, we consider the mappings <J> and ^ = 
defined by 

*(Q)(x) = Q[^(x)] and $(P)(u) = P[^(u)] , 
where P £ V (A n>r ) and Q € V (Pv,„). 

Remark 2.1. Note that P is the uniform distribution on ^4„ jr (i.e., the Bose-Einstein model) 
if and only if Q = <3?(P) is the uniform distribution on Pv 5 „. 

Consider now r exchangeable random variables Y±, Y2, . . . , Y r , which take values in {1, 2, . . . , n}. 
To these random variables we can associate a vector of occupancy numbers by introducing 
the random variables X± , X2 , . . ■ , X n defined by 

r 

X 3=J2 1 {Y h =j}, for j = 1,2,..., n. 
h=i 

The set of the possible values taken by Y = (Y\, . . . , Y r ) is then D r ^ n := {1, 2, n} r . 
We will also write X = <p (Y) or 

X j = <p j (Y 1 ,Y 2 ,...,Y r ), forj = l,...,n (2) 

where (p : D rj7l — > A n ^ r with 

r 

<Pj{yi,V2, ■ ■ ■ ,Vr) = ^Z 1 {y h =j}i for y G I? r , n and j = 1,2, ...,n. (3) 
h=i 

It can be immediately seen that the probability distribution of (Yi, Y2, . . . , Y r ) is uniquely 
determined by the probability distribution of (X\,X 2 , . ■ ■ , X n ) and vice versa. 
In fact, the following relationships hold 

vrv v v 1 P{X 1 = <p 1 (y),X 2 = (p 2 (y),---,X n = (p n (y)} 

P{Y\ = yi, Y 2 = y 2 , ■ ■ ■ , Y r = y r } = r (4) 

Wi(y)^2(y)--^n(y)/ 

and 

P{Xi = xi,X 2 = x 2 , ■ ■ ■ ,X n = x n } = 

r )p{y 1 = ^(x),y 2 = v 2 (x),...,i; = Vr(x)} (5) 

X\ • • • X n J 

We note that (p : D r ^ n — > A n , r is different from the previous transformation ip, which is 
defined on B r ^ n C T) r ,n an d is bijective. However, <f(y) = y(y) if y G B r ^ n and, for any 
y G D r ^ n , we have 

= <p(V(i), %)) = <p(V(i), V(r)), 
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where (2/(1), • • • , Um) is the vector of the coordinates of y rearranged in increasing order. 

We consider now the vector Yr,\ = (Yn\, . . • , Y( r )) °f the order statistics of the exchangeable 
vector (Yi,Y%, . . . ,Y r ). The set of values taken by Y(.) is B r>n . To a probability distribution 
of (Yi,Y2, ■ ■ ■ ,Y r ) it corresponds one and only one probability distribution of Yq and, for 
u € B rtn , we can write 

\ipi(u)—tp„(u)) 

Furthermore, 

<^(Y ( . ) )=^(Y ( . ) )=^(Y) = X, 

and thus 

Y ( . ) =(^- 1 (X)=V(X). 
We can then summarize the arguments above as follows. 

Proposition 2.1. Let X\, X2, ■ ■ ■ , X n be random variables such that the support of their 
joint distribution is A n ^ r , for some r € N. Moreover, let Yi, Y2, ■ ■ ■ , Y r be exchangeable, 
{1, 2, . . . , n}-valued random variables such that ^ holds. 

Then, the random variables Ui, U2, ■ ■ ■ , U r , defined by equation ([T]), are the order statistics of 
the vector (Y\ , . . . , Y r ) . 

In what follows we present an interpretation of the marginal distributions of the variables 
Y%, Y2, . . . , 1^., in terms of the associated occupancy numbers X\, X2, ■ ■ ■ , X n . On this purpose, 
we now consider a special transformation, that map one occupancy model into another. Two 
other natural transformations, also intersting for our next developments, will be introduced 
later on. 

2.1 Transformation JCi 

The first transformation we consider can be described as follows. 

Consider r objects distributed on n cells according to a given occupancy model and let 
Xi,X2,.. . ,X n be the related occupancy numbers. We drop one of the particles from this 
population randomly (i.e. so that each of the r particles has the same probability to be 
dropped, independently of its cell). We denote by X[,X' 2 , ■ ■ ■ ,X' n the occupancy numbers 
associated with the new population of r — 1 particles. 

This transformation will be denoted by fC\. Notice that K,\ : V (A nr ) — > V (A nr —i). 
Lemma 2.1. Consider the event E^i = {X[ = x\, . . . , X' n = x' n }. Then, 

P{E X ,} = £ X ^P{X l =x' 1 ,...,X h = x' h + l,...,X n = x'J. (6) 

h=l r 
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Proof. The proof consists of simple manipulations. In fact, we can write 
P{E*} = P{X[ = x' 1 ,...,X' n = x' n } 

n 

= P{E X > H (the eliminated particle is from the cell h)} 
h=i 

n 

= Y J P{E* n (X 1 = x' 1 ,...,X h = x' h + l,...,X n = x' n )} 

h=l 

n 

= Y J P{E X '\X 1 = x' 1 ,...,X h = x' h + l,...,X n = x'Jx 
h=i 

x P{X! = x[,...,X h = x' h + l,...,X n = x' n } 
n ' + 1 

= — P ^ Xl = x i> ■■■,Xh = x' h + l,...,X n = x' n }. 



r 
h=l 



□ 

Let us now consider the exchangeable vectors (Y\, . . . , Y r ) and (Y{, . . . , Y/_]J , corresponding 
to the occupancy numbers X\, . . . , X n and X[, . . . , X' n , respectively. 

Proposition 2.2. The joint distribution of (Y[, . . . coincides with the (r — 1)- dimen- 

sional marginal distribution of (Y\, . . . , Y r ). 

Proof. In view of Q and Q, the joint distribution of (Y/, . . . , Y 7 f_ 1 ) is given by 

P { Y { = Vl, Y 2 = 2/2> • • • , Y r-1 = y'r-l) = 



P{X[ = gifrQ, X> 2 = y 2 (y'), . . . , = 
Ui(y')*a(y , )-^»(y')) 



= yi(yQy2(yQ • • • yn(y') £h(y') + 1 

(r-l)! ^ r 

x P{X l = My), ■■■,X h = <p h {y') + l,...,X n = <p n (y')}. 

For brevity sake, let us set 

^W(y') := (v3 1 (y0,...,^(y0 + l,...^n(y')) , 
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so that we rewrite the previous equation as 



P{Y{ = y[,Y£ = y' 2 ,..., = y' r _ x } = 

= ^i(y / )^2(y / )---^n(y / ) A >Ph(y') + 1 P r x = ^OOfv'n 
(r-l)! 2^ r i ^ ^ y ^' 

We can now use ^ and get 



^i(y')!^(y / )!--^n(y , )!^^(y / ) + i H 



(r-i)! ^ r ^i(y')!- ••fe(y') + !)!■•• ^n(y')! 
x p{y x = (V h) (y')) = Vv-i (^(y')) ,n = Vv (V^V))} 

n 

= E P = ^ {^ h) (y')) > • • • , = W-l (<f (h) (y)) , Y r = t/> r (<P (h) (/)) } • 

Let us focus on the vector 

(v^V)) ,Vv (V h) (y'))) • 

The first r — 1 components are defined by the composition (-0 o (p){y') and the result is a 
(r — l)-tuple 

(y'a(l)' ■ ■ ■ >y'a{r- 1)) € ^r-l.n 

for every permutation cr of {1, 2, . . . , r — 1}. The last component ip r (^^ h '(y')) is n °t deter- 
mined by the previous composition. It is a "free" element indicating the cell in which the 
additional r-th particle fell. Thus, it is a value randomly chosen in {1, . . . , n}. 
Finally, by recalling the exchangeability of (Yi, . . . , Y r ), we get the equality 

n 

P{Y{ = = y' 2 , ■ ■ ■ , Y' T -X = V'r-l} = = y'l' ■ ■ • ' ^-1 = y'r-l,Yr = h}. 

h=l 

and then the conclusion follows. □ 



2.2 Transformation /C 2 

Here we consider the transformation K-2 '■ V (^4 n ,r) — ► P (-<4n-i,r) simply obtained by "eras- 
ing" one of the n cells. More precisely, we start once again from an occupancy model in 
V (A njr ) and let X\, . . . , X n denote occupancy numbers jointly distributed according to such 
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a model. We now consider the case where the n-th cell is eliminated and any of the X n 
particles that had fallen within it, is put at random, and independently, within the remaining 
cells (see also details about the MB model later on). 

By denoting X^,...,X^_ 1 the occupancy numbers obtained by this procedure, and for 
(x'(, . . . , x'n-i) € A n -i !r , we can write 

P{X'( = x", . . . , X'n_! = a;"_i} = 

( X ) 

= E E (n-lY P{Xl = X "-^---> X n~l = <-l - Xn = x}, (7) 

where we obviously mean P{X± = x±, . . . ,X n = x n } = 0, whenever some of the coordinates 
x±, . . . , x n _i turns out to be smaller than 0. 

2.3 Transformation /Ci^ 

Here and in the rest of the paper we will use the notation SV = J2iLi We consider the 
transformation that maps V (An :J .) onto V (A njS ) (with 1 < n < N — 1,1 < s < r) that 
is simply obtained by conditioning on a fixed value s for the partial sum S n . For variables 
Xi, . . . ,Xn distributed according to a given occupancy model in V (A/v,r) we consider the 
model in A n ^ s defined by 

vfv Y \a \ P{Xi=x 1 ,...,X n = x n } 
P{Xi = xi, . . . , X n = x n \S n = s\ = _ 

*\&n — s f 

for (xi, . . . , x n ) G A n<s and with 

P{X 1 =xi,...,X n = x n } = 

= ^ P{X 1 = x 1 ,...,X n = x n ,X n+ i = r] 1 ,...,X N = r) N - n }, 

— n,r — s 

P{S n = s}= ^ ^ P{X 1 =X 1 ,...,X n = X n ,X n+1 =7] 1 ,...,X n = 7] N . n }. 

— n,r — s 

3 Exchangeable Occupancy Models 

Often relevant occupancy models are such that the random variables X\ , . . . , X n are ex- 
changeable. The class of the Exchangeable Occupancy Models (EOM) is actually a wide and 
interesting one. We then devote this Section to analyze several special aspects related with 
such a condition. We notice, in particular, that the most well-known occupancy models, i.e. 
Maxwell- Boltzmann, Bose-Einstein and Fermi-Dirac, are exchangeable. They are defined as 
follows. 
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• Maxwell-Boltzmann: 



P{Xi = Xl,X 2 =X2,-..,X n = x n } 



n r xi\x 2 \ • • • x n \ ' 



for x = (xi,x 2 , . . . , x n ) G A n , 



Bose-Einstein: 



P{Xx = xi,X 2 = x 2 ,...,X n = x n } = 

{ n-l ) 



for x = (xi, x 2 , ■ ■ ■ , x n ) G A n . 



Fermi-Dirac: 



P{Xi = xi,X 2 = x 2 , . . . ,X n = x n } = -py 



for x = (x\,x 2 , . . . , x n ) G A n , ri where A n>r := {x E A n>r : Xj G {0, 1}}. 

As mentioned above, we can easily see that any vector (X\, . . . , X n ) distributed according to 
one of these models is an exchangeable random vector. 

In the next Section we will introduce a wide class of EOM's that contains these fundamental 



models. For an account about MB, BE, FD see Feller (1968) and Charalambides (2005). In 



the following Propositions we state some general properties of EOM's. 

Proposition 3.1. If X\, X 2 , . . . , X n are exchangeable, then the univariate marginals ofY\, 
Y 2 , . . . , Y r are uniform on {1,2, ... ,n}. 

Proof We prove the statement by showing that for any a, b G {1,2,..., n}, a ^ b, we have 
P{Y\ = a} = P{Y\ = b}. This would mean that P{Y\ = a} = ^ and hence the univariate 
marginal distributions of (Y±,Y 2 , . . . , Y r ) are uniform on {1, 2, . . . , n}. 
Let (y 2 , . . . , y r ) be a (r — l)-tuple in D r _i jn and consider a G {1, 2, ... , n}. Then 

P{Y 1= a}= P{Y 1 =a,Y 2 =y 2 ,...,Y r =y r }. 

(j/2,---,J/r)eA'-l,n 
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By using Q we obtain that 

P{X 1 = <pi(a,y 2 , ■ ■ .,y r ), ...,X n = <p n (a,y 2 , ■ ■ .,2/r)} 



PiY, = a}= Yl 



(lAi,— ilfcOeA— l,n \-Vl(0)10,»-,J/r)^2(o,J/2,...,»r)— <fin(a,y2,-,Vr) 

<pi(a, y 2 , ■ ■ ■ , 2/r)! • • • <p n (a, y 2 , ■ ■ ■ , VrV- 



E 



■X 



r! 

(B,-,!/r)6fr-l,n 



x P{Xi = (p\{a,y 2 , . . .,y r ), ...,X n = (p n {a,y 2 , . . . ,y r )} 

Taking into account that for any j G {1, 2, . . . , n}, the function ipj is given by ([3]), we obtain 

Dfv v- tpx(y 2 ,..., y r )\ ■ ■ ■ (<p a (y 2 , . . . ,y r ) + 1)!-- -<p n (y 2 ,..., y r )\ 
P{Yx = a}= ^ x 

(l/2,...,!/r)eDr-l,n 

x P{X X = <pi(y 2 , . . .,y r ), ...,X a = <p a (y 2 , • • • , y r ) + 1, ■ ■ ■ , X n = tp n (y 2 , • • -,2/r)} 

E<fi(y2,---,yrV----<fn(y2,---,yr)K~ , . , ^ 
-j (<rM2/2, • • ■ ,2/r) + 1) 

fc,.~,!/r)EOr-l,» 

x = i^i 0/2, • • • ,2/r), ...,X a = (p a (y 2 , . . . ,y r ) + I, ■ ■ ■ , X n = (p n {y2, • • -,2/r)}- 

(8) 

Now, we take 6 G {1,2,..., n}, b ^ a. Without loss of generality we can suppose that a < b. 
Then, the univariate marginal of Y\ computed in b is 



P{Yl=6}= 2^ (mi/2,...,i/r) + i) 

(y2,-*)6-D r _i,„ 

X P{Xi = ^1(2/2 , • • .,2/r)) • • • = *Pb(y2, ■ ■ -,2/r) + 1, • • ■ ,X n = <p n (y2, ■ ■ -,2/r)}- 

(9) 

The sums ([8]) and ([9]) have n r_1 terms and are computed using all possible (r — l)-tuples 
(2/2 j • • • j 2/r) £ P > r-i,n- Therefore for any j G {1, 2, . . . , n} and i G {0, 1, 2, . . . , r — 1} there are 
summation terms in ^ and in Q, such that <Pj(y 2 , . . . ,y r ) = i. Let a, (3 G {0, 1, . . . , r — 1} 
two fixed values. Then there are summation terms in ^ such that (p a {y2, ■ ■ ■ , 2/r) = a and 
<?&(2/2, - - • , 2/r) = /3- Also there are summation terms in ([9]) such that (p a (y 2 , . . . ,y r ) = (5 and 
<^ft(2/2, ... ,2/r) = a- 

Let us now consider (xi, . . . , x a , . . . , x b , . . . , x n ) G -A n ,r-i such that x a = a and x b = /3. Then, 
there are f £~ „ ) summation terms in (f8l) , all equal to 

]_ ^ * * * ! 

— ' — : ' (a + l)P{Xi = x 1} . . . ,X a = a + 1, . . . ,X b = (3, . . . ,X n = x n } 

r\ 
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The overall sum of these terms is: 
a + 1 



P{Xi = xx, . . . ,X a = a + 1, . . . ,X b = P, . . . ,X n = x n }. 



(10) 



In the same way let us consider (x±, . . . , x a , . . . , Xf,, . . . , x n ) G A njT ._i such that x a = j3 and 
Xb = a. Then, there are f iT 1 ^. ) summation terms in (181), all equal to 



(a + l)P{Xi = a?i, . . . , X a = P, ■ ■ ■ ,X b = a + 1, . . . , X n = x n } 



The overall sum of these terms is: 
a + 1 



P{X X = xi, . . . ,X a = P, . . . ,X h = a + 1, . . . ,X n = x n }. 



(11) 



The random vector (X\,X2, ■ . ■ ,X n ) is exchangeable and therefore (10) and (11) are equal. 



For any a, P G {0, 1, . . . ,r — 1} the above terms uniquely determine the sums © and ©, 
which then are equal. As a consequence, P{Y\ = a} = P{Y\ = b}. 

□ 



Remark 3.1. Let 



(n) r = n(n + 1) ■ • • (n + r — 1) 



n 



(r) 



n(n — 1) • ■ ■ (n — r + 1) 



be the ascending and falling factorial, respectively (see Charalambides (2005)). For MB, BE 



and FD, the joint distributions of the random variables Y±, . . . Y r are given by 
• Maxwell-Boltzmann: 



P{Y 1 =y 1 ,...,Y r = y r } 



1 



TV 



for y G D r n ; 



Bose-Einstein: 

P{Y 1 =y 1 ,...,Y r = y r } 
for y G D r n ; 



<^i(y)!---^n(y)! 



(n + r — l)(n + r — 2) • • • (n + l)n 



n?=ifr(y)' 

(n) r 



Fermi-Dirac: 



p{y x = yi,...,y r = y r } 



<^i(y)!---0n(y)! 



n-=i^(y)! 



n(n — l)(n — 2) • • • (n — r + 1) 



n 



(r) 
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for y G D r>n , where D rjU = jy G -D rjn : X)j=i — n an d the y/s are all distinct j. 



All these three distributions admit uniform marginals, in agreement with Proposition 3.1 



Observe that, in the MB case, the random variables Yi,...,Y r are independent and their 
distribution is uniform over the set {1, 2, . . . , n}. 

As a remarkable property of the class of the EOM's we have the closure under the transfor- 
mations Ki, K.2 and K^f\ 

Proposition 3.2. Let (X\, . . . , -Xjv) be an EOM on A^^ r and fix n < N. Conditionally on 
the event {S n = s}, the variables X\, ■ ■ ■ ,X n are distributed according to an EOM on A n ^ s . 

Proof. It must be proved that conditionally on {S n = s}, the vector {X\, . . . , X n ) is exchange- 
able, i.e. 

P {^(l) = xi,.. .,X a{n ) = x n \S n = s) = P{Xi =xi,...,X n = x n \S n = s} 

for every permutation a of {1, 2, . . . , n}. It straightforwardly follows from the exchangeability 
of the random variables X±, . . . , X^. Indeed, for any {x\, . . . , x n ) G A n<s , we obtain 

=xi,.. .,X a ^ = x n \S n = s} = 

_ P{X\ = X a [i) , X a ( n ) = x n } 
P{S n = S} 



E 

— n,r — s 



P{Xi = x a (y, . . . , X a ( n ) = x n , X n+ i = rjx, . . . , X N = n N - n } 

P{S n = s} 



E 

— n,r — s 



P{Xi = xi, . . . ,X n = x n ,X n+ i = rji, . . . ,X N = r]N-n} 
P{S n = s} 



P{X 1 = x 1 ,...,X n = x n , S n = s} 

P{S n = S} 



= P{Xi = xi, . . . , X n = x n \S n = s} 
and this concludes the proof. 



□ 



The following Proposition can be easily derived from Lemma 2.1 
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Proposition 3.3. Let (X\, . . . , X n ) be an EOM on A nr and (X[, . . . , X' n ) be the occupancy 
numbers obtained by applying the transformation 1C\. Then, (X[, . . . , X' n ) is an EOM on 

-^n,r—l ■ 

The following Proposition can be easily obtained by taking into account equation Q. 

Proposition 3.4. Let (X\, . . . , X n ) be an EOM on A n ^ r and (X[, . . . , X' n ) be the occupancy 
numbers obtained by applying the transformation K,<i- Then, (X[, . . . , X' n ) is an EOM on 
A-n— i,r • 

Remark 3.2. The transformation K2 can be seen as a special case of a more general class 
of transformations: we can consider the case where the n-th cell is eliminated and the X n 
particles that had fallen within it are distributed within the remaining cells according to an 
exchangeable occupancy model. The class of the EOM's is closed also w.r.t. this class of 
transformations . 



4 Occupancy .M ^-models 



In this Section we consider a remarkable sub-class of EOM's, that includes the models MB, 
FD and BE. Such a class can be introduced as follows: fix a function a : {0, 1, . . . } — > R+ 
and for r, n G N, set 

n?-i a(xj) 

P{X = x} = J (a) for x G A n , r , (12) 

where 

n 

£eA n>r j=i 

In this case we say that X = (X\, . . . , X n ) is distributed according to the exchangeable occu- 
pancy model Jv6n)r- 

Notice that the Maxwell-Boltzmann, Bose-Einstein and Fermi-Dirac models are respectively 
obtained by letting 

MB: a(x) = -r, for x = 0, 1, 2, . . . 
x\ 

BE: a(x) = l, for x = 0, 1, 2, . . . 

FD: o(0) = 1, o(l) = 1, a(x) = 0, for x = 2,3,... 



It is interesting at this stage to point out the existence of .M^r-models different from the three 



above. An example is provided by the Pseudo-contagious models presented in |Charalambides 

npoB). 
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Example (Pseudo-contagious occupancy model). This is the model characterized by the joint 
probability distribution 

^S+Xl-l^ _ _ _ fs+x„ — l\ 

P{Xi =xi,...,X n = x n } = Xl , sn+r _ u Xn , for x G A n , r . 



^sn+i — 1^ 



This model belongs to the class .M^r- In fact, it corresponds to the choice of a(x) = ( s+x . 1 )- 

In a A^r-model, the joint distribution of the Y\, . . . , Y r random variables introduced Section 
[His 

P\ Y 1 — Vl, y 2 — 2/2, • • • , Yr — Vr) — 73 • 



Charalambides 



(2005, Chapter 4) 



An approach for constructing A^r-models is described in 
and makes use of i.i.d. random variables. It works as follows. Let Z\, . . . ,Z n be N- valued 
i.i.d. random variables with common law 



P{Zi = x} = q(x) for i = 1, . . . , n. 



and let S n = Z\ H + Z n . 

Consider the random vector (X±, . . . , X n ) with values in A n%r and such that 

P{X 1 =x 1 ,...,X n = x n } := P{Z 1 =xi,...,Z n = x n \S n = r}. 
This means that, for any (x±, . . . , x n ) 6 A n>r , it holds 

P{Z 1 = Xl, . . . ,Z n = X n } 



P{Xi = Xl, . . . ,X n = X n ] 



P{S n = r} 



P{S n = r} ' 



(13) 



The Formula (12) can be seen as a slight generalization of Formula (13); actually a(x) is not 



necessarily a probability distribution. We can furthermore extend as follows the construction 
of the approach based on i.i.d. variables. Let Z\,...,Z n be N- valued conditionally i.i.d. 
random variables with joint discrete density 



" r+ao 

P{Zi = zi, . . . , Z n = z n } = Y[ a (zj) / exp < 

j=i J -°° 



-Oj^zj }A(d6) 

3=1 



(14) 



A being a probability distribution on the real line. Put S n = Y^j=i %j an d let us repeat the 
same construction as before: consider the random vector (X\, . . . ,X n ) with values in A n ^ r 
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and such that, for x S A nr , 

P{Xi = xi, . . . ,X n = x n } := P{Z\ = xi,. . . ,Z n = x n \S n = r}. 



Then, as it is easily seen, we reobtain the occupancy model in (12). Indeed, we have 



n /7 ry \Q \ P {^1 = Xl, . . . , Z n = X n ] 

P{Z X =xi,...,Z n = x n \S n = r\ = — — 

P{S n = r] 

n" = i"(^)L + ^exp{-gr}A(^) 

E ve ,„, r n u a w *M-er}A m 
n"=i«(^) 

r<( a ) 

Remark 4.1. The fact that P{Z\ = x%, . . . , Z n = x n \S n = r} does not depend on the distri- 
bution A, has an immediate interpretation in statistical terms: for the exponential model in 



(14), S n is a sufficient statistic with respect to the parameter 9. 



The construction considered above is somehow more general than the one based on i.i.d. 



variables. Anyway, the Formula (13) suggests some important remarks. 



Remark 4.2. Up to a multiplicative factor, the normalization constant C^l in (12) can be 
interpreted as the probability that a sum of certain i.i.d. random variables is equal to r, i.e. 
Cn"r = kP{S n = r} with k suitable constant. 

Remark 4.3. With the above construction the three fundamental models can be recovered 
by choosing suitably the distribution q(x). In fact, MB, BE and FD correspond to the cases 
when q(x) obeys a Poisson, Geometric and Bernoulli distribution, respectively. Moreover, the 
pseudo-contagious occupancy model is obtained by setting q(x) a negative Binomial distribu- 
tion. 



From Proposition 3.2 we know that a subvector (Xi, . . . , X n ) of an EOM (Xi, . . . , Xn) is still 



an EOM conditionally on a fixed value for the sum S n (closure property w.r.t. transformations 
of the type K,^/^). Then we know that, by applying this transformation to an occupancy 
model of the form A4^ a \ we surely obtain an exchangeable model. It is then natural to wonder 
whether the latter is still of the form Ai^ (i.e., whether the closure property also holds for 
the class Ai^). The following proposition answers this question. 

Proposition 4.1. Let {X\, . . . , Xn) be a M$ -model and define S n = X\ + • • • + X n with 
n < N. Conditionally on the event {S n = s}, the variables X\, ■ ■ ■ ,X n are distributed as a 
M-ns-fnodel. 
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Proof. We need to prove that conditionally on {S n = s}, the vector (X\, . . . , X n ) is distributed 
as a A^n^-model, i.e. 



nfy V | 1 Ylj = i a ( X j) 

P{Xi = xi, . . . ,X n = x n \S n = s\ = — 7-t 



For any (xi, . . . ,x n ) £ A n ^ s , we have 

P{Xx = x\, . . . ,X n = x n \S n = s} = 

P{Xi =xi,...,X n = x n } 



E 



P{S n = S} 

P{X 1 = X 1 ,...,X n = X n , X n+1 = 771, . . . ,X N = T] N -n} 



P{S n = S} 

^ U]=i a ( x j)U^=i n aim) 

r?eA JV _ nir _ s C^ r P{S n = s} 

n 

KY[a( Xj ) 

3=1 



where 



N-n 



k ■= ^) „ r „ , e n «*) ■ 



C Ntr P{S n - s} ve A N _„ tr _ 3 i=l 

Since P{X\ = x%, . . . ,X n = x n \S n = s} is a probability distribution, K is exactly equal to 
{Cnls^j and this concludes the proof. □ 



As seen in Proposition 3.3 the action of dropping one particle at random (transformation 
tCi) preserves the exchangeability property. As we will see in next Proposition, the structure 
of A4^ -model is preserved only under the technical assumption (15). It is easy to see that 



MB, BE, FD and the pseudo-contagious occupancy models satisfy this condition. 

Proposition 4.2. Let (Xi, X%, . . . , X n ) be a Aify-model and let the function a : {0, 1, . . . } 

satisfies the condition 



C n)--i x' h + 1 a(x' h + 1) , 

^n,r h=l y h > 

Then, the occupancy vector (X[, . . . ,X' n ), obtained by applying the transformation fC\, is 
distributed according to the model 
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Proof. We have to show that 

nj-i a K 



P{X[ — x 1 , . . . , X n — x n } 



u u,r-l 



In view of ^ the joint distribution of (X[, . . . , X' n ) is given by 

n x' + 1 

= x[, . . . , X' n = x' n } = ^ — P{X X =x' 1 ,...,X h = x' h + l,...,X n = x' n } 



h=l 

h=l ^n,r j=l 

n^ 1 «K)gg- 1 y. < + laK + l) 
^ r a(x' h ) 



where the last equality follows from (15). □ 



Remark 4.4. One can easily realize that the class of the Ai^ -models is strictly contained 
within the class of the EOM's. Examples can be easily found, for instance, by starting from 
models in the class M.^ and by applying the transformations K,\ or K-2- 

5 Processes with the generalized UOSP 

For M G NU{+oo}, let {Nt}t=o,l,...,M be a discrete-time counting process with jump amounts 
Jo, Ji, . . . , Jm, i-e. Jo, Ji, • • • , Jm is a sequence of {0, 1,2,... }-valued random variables and, 
for t = 0, 1,...,M, 

t 

N t = J2 J h- ( 16 ) 
h=0 

For our purposes we introduce the following notation and definitions. 
Let Ji,T2, . . . denote the arrival times of the process {N t }t=o,i,..., i.e. 

T n = mf{t >0:N t >n}, (17) 

and let Zx, Z2, ■ ■ ■ denote the inter-arrival times, i.e. 

Zn = Tn Tfi—X • 



17 



Notice that, since P{Jh > 1} > 0, it can happen {T n ^\ = T n }, and then {Z n = 0}, for some 
n. The zeros of the Z's are related to the ties of the T's and to the jump amounts greater 
than one. 

Definition 5.1. Let {Nt}t=o,i,... be a discrete-time counting process with jump amounts 
Jo, ... , Jt- We say that it satisfies the -Uniform Order Statistics Property (.A/f ^-UOSP) 
if, for any t, k G N and any (jo, ■ ■ ■ ,jt) £ ^t+i,k, we have 



P{J =j ,...,J t =j t \N t = k} 



nLo a (^) 



(a) 

t+l,k 



(18) 



This definition can be seen as a natural, and unifying extension of the definition of discrete 



UOSP given in |Shaked, Spizzichino and SuterJ ( [2004J [2008J ) . 

Definition 5.2. Let a : {0, 1, . . . } — > M + and Rt : {0, 1, ... } — > M + be given functions. 
The process {N t }t=o,i,...,M is an a-mixed-geometric process if the discrete joint density of 
(Jo, Ji, . . . , Jt) has the form 



Pt(jo,jl, ■ ■ ■ , Jt) = Rt[ z^jh ■ a (^) for t = 0, 1, . . . ,M. 



(19) 



\h=0 



h=0 



Remark 5.1. The sequence of functions R\, R2, ... in the Definition 5.2 cannot be independent 
of function a. In fact, since the discrete density pt-i must be the marginal of pt, from (19) 
we obtain 



Rt-i(k)=J2a(l)R t (k + l). 



1=0 



We can now formulate the following Theorem, that extends the results in Shaked, Spizzichino 



and Suter (2004 Section 4). 



Theorem 5.1. Let {iVt}t = o,i,... be a discrete-time counting process. Then, the following 
statements are equivalent: 

(i) {iV t } t= o,i,... satisfies the M^ a) -UOSP. 

(ii) {Nt}t=o,i,... is an a-mixed-geometric process. 

(Hi) For any k E N, the discrete joint density of the inter-arrival times Zi, Z2, ■ ■ ■ , Zf. is 
given by 

£<=i «. / \ 
P{Z 1 = z 1 ,...,Z k = z k } = R EUzi (k)- J] ffi E^Ek^} ' ( 2 °) 

h=0 \ i J 



18 



(iv) For any k G N, the discrete joint density of the arrival times T%, T2, ■ ■ ■ , T\~ is given by 

tk / \ 



P{T X = tl ,...,T k = t k } = R tk (k) • II ° E 1 W= 



h} 



(21) 



h=0 \ i 



Proof, (i) => (ii). Consider {N t } t =o,i..., satisfying the .M^-UOSP. For any k G N, we set 

P{N t = k} 



Rt(k) r- 



a 



(a) 

t+l,k 



Clearly (19) is satisfied. In fact, it easily yields 



P{Jo =jo,---,Jt = jt} = P I Jo = jo, ■ ■ ■ , Jt-i = jt-i,N t = ^jh 

I h=o J 

iVt = X>| p\N t = Y J 3h 

h=0 ) I h=0 > 



^ j^O = jo, h- 1 - .//-i 

nLo a 0>o 



c 



(a) 



fc=0 



where the last equality is due to (18). 



(ii) (i). Consider {Nt}t=o,i..., an a-mixed-geometric process. For any t, k E N and for any 
(jo, • • • Jt) € ^t+l,fc, by ([T9| we get 

^{ Jo = jo, ■ ■ • , Jt = jt} 
P{N t = k} 

M k ) ■u t h=o a Uh) = nLo a Uh) 



p{j =io,...,Jt = i*|JVt = fc} 



which is exactly (18). 



(ii) =>• (iii). Consider {Nt}t=o,i..., an a-mixed-geometric process and for any k G N, let 
Z\,...,Zk be positive integers. We have 



P{Zi = zx, . . . , Z k = Zk} = P <Ti = Zi,T 2 = zi 



k 1 



Moreover, by the definitions of counting process (16) and of arrival times (17), we obtain that 
the jump amounts are given by 
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Hence, we get 

P j^i = zi, T 2 = zi + z 2 , . . . , Tfc = ^2 2 d j = 

= p \ J o = Yl 1 m=o},Jz 1 = %=«i}' • • • ' J eLi «a = zZ V*=ELi f 

I i i i ) 

Ed=i *d / \ 
= ^Lx*( Jo + ,7 * 1 + *" + J £5-x*J II a E^J 

h=Q \ i / 
E5=i «i / \ 

= /2 EJU* i W II ° I «*=*}) 

h=0 \ i / 

and the implication is proved. 

(iii) =>• (ii). Let L = Ya=o l{Ji>o}! i- e - the number of J's that are strictly positive. Define 

h := inf{0 < i < t\J t > 0} 

I n : = inf{/ n _i < i < t\Ji > 0} , 1 < n < L 

and let 

n 

J n :=^J/ h n = l,2, ...,L. 
h=l 

Note that the sequence I\, . . . ,Il is strictly increasing and moreover, 



L t 

>b ■ 



Jl = Y Jl h = E Jh 



h=l h=0 

Without loss of generality, suppose that { Jt > 0} . Hence II = t and we have 

Z 1 =h Z 2 =0 ••• z 3l =0 

z Jl+1 =h-h z Ji+2 =0 ••• Zj 2 = 



^Ji-l+l =t-I L -l Z J L _ 1+ 2 =° Z J L =°- 



(22) 



Prom (22), the joint distribution of (Jo, Ji, . . . , Jt) can be written as 



P{Jo = jo, Ji = 31, ■ ■ • , Jt = jt} = P{Zi = h, Z 2 = ■ ■ ■ = Zj x = 0, Zj i+1 =i 2 -h,... 



Z J L - 1 +i = 1 ~ 'i-i'^+J = --- = Z J L =°} 
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Since 

ji 

I> = * and E ^Eki **=h} = Oh 



i=l 



from (20) it follows 

P{Jo = Jo, h =ji,---,Jt = jt} = Rt I ^2,3h I J| a (j h ) . 

\h=0 / h=0 

(iii) => (iv). For any k £ N, consider ti, . . . , positive integers. We can write 
P{Ti = ti, . . . , T k = t k } = P{Z\ = ti, Z2 = t2 — ti, . . . , Zj, = tfc — ife-i} 

= R T,Li ^-td-i n ° e i {Ei=i t d -*d-i=h} 



11° E x m= 



h=0 



(iv) =^ (iii). For any £ N, consider Zi, . . . , z k positive integers. We have 

r k 

P{Z X = zi, . . . , Z k = z k } = P I Ti = zi, T 2 = zi + z 2 , ■ . . , T k = z d 



d=l 



and the conclusion immediately follows by (21) 



□ 

Here, we conclude the paper with a few comments about the processes with the .M^-UOSP. 
As it is clear from the definition, the role of such a property is played around conditioning 
w.r.t. events of the form {Nt = k}. 

Let us consider, in the beginning, a general sequence of exchangeable jump amounts J\,J 2 ,... 
and put N s = Y2t=i Ji, s = 1,2, ... . Then, conditionally on {N = k}, the joint distribution of 
J\, J2, ■ ■ ■ , Jt is an EOM over A tj k- Let Y-[ , . . . , 5v be the {1,2, ... ,t}- valued, exchangeable, 
random variables that correspond to such an occupancy model. We denote by their joint 
probability law and recall that their common marginal distribution is, in any case, uniform 
over {1, 2, . . . ,t}. 

In the special case of processes {N s } s=h2 ,... with the A^W-UOSP, £ (t ' fc) is given by 

p \y} t] = yi, y® = y 2 ,..., if = Vk } = nL oQ (^(y))^(y)! . 



kicr + \ 
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When more in particular a{x) = \/x\, i.e. when the conditional distribution of J\, J2, ■ ■ ■ , Jt 
is the MB model over A t ^, . . . , are i.i.d. with uniform distribution over {1, 2, . . . , t}. 
This case, the one denoted by UOSP(<) in Shaked, Spizzichino and Suter ( 2004[ ), can be 
seen in a sense as the discrete-time analog of the homogeneous Poisson process (in continuous 
time). The latter satisfies in fact the classical Order Statistics Property: conditionally on 
{Nt = k}, the arrival times T%,. . . can be seen as the order statistics of k independent 
random variables, with an uniform distribution over the interval [0, t]. The set of the arrivals 
can thus be seen as a random sample from a uniform distribution over [0, t]. In the case of 
UOSP(<) then we have that, for any k > 1 and t > 1, obviously coincides with the 

(k — l)-dimensional marginal of £^' k \ This last circumstance actually is a weaker condition 
than UOSP(<), in fact it holds for any .M^-UOSP process with a satisfying the condition 
(15). It can however be still considered as a condition of randomness for the arrivals, in view 



of the meaning of the transformation /Ci and of Proposition 2.2 
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